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ABSTRACT 

Measuring the distribution of mass on galaxy cluster scales is a crucial test of the 
ACDM model, providing constraints on the nature of dark matter. Recent work in- 
vestigating mass distributions of individual galaxy clusters using gravitational lensing 
has illuminated potential inconsistencies between the predictions of structure forma- 
tion models relating halo mass to concentration and those relationships as measured 
in massive clusters. However, such analyses typically employ only simple spherical 
halo models with canonical NFW slopes, while the halos formed in simulations show 
a range of more complex features. Here we investigate the impact of such expected 
deviations from the canonical NFW halo profile on mass and parameter estimation 
using weak gravitational lensing on massive cluster scales. Changes from the canoni- 
cal NFW profile slopes are found to affect parameter estimation. However, the most 
important deviation is halo triaxiality because it is impossible even with fiducial weak 
lensing data to fully resolve the three-dimensional structure of the halo due to lens- 
ing's sensitivity only to projected mass. Significant elongation of the halo along the 
line of sight can cause the mass and concentration to be overestimated by as much as 
50% and by a factor of 2, respectively, while foreshortening has the opposite effect. 
Additionally, triaxial halos in certain orientations arc much better lenses than their 
spherical counterparts of the same mass, indicating that clusters chosen for study be- 
cause of evident lensing are likely to be drawn from the high-triaxiality end of the halo 
shape distribution; cluster samples chosen with no shear bias return correct average 
parameter values. While the effects of triaxiality alone may not be enough to fully 
explain the very high concentrations reported for some clusters, such as Abell 1689, 
they go a long way in easing the tensions between observations and the predictions of 
the cold dark matter paradigm. 

Key words: gravitational lensing - cosmology: theory - dark matter - galax- 
ies: clusters: general. 



1 INTRODUCTION 

Galaxy clusters are ideal laboratories in which to study 
dark matter, being the most massive bound structures in 
the universe and dominated by their dark matter compo- 
nent (~ 90%). Constraining the clustering properties of 
dark matter is crucial for refining structure formation mod- 
els that predict both t he shapes of dark matter halos and 
their mass function (e .g. iNavarro" et al.l(|l997l '): lBahcall et al.l 
l|2003h : iDahld (|2007l )). Several methods are used to mea- 
sure galaxy cluster dark matter profile shapes and halo 
masses on a range of scales, including X-ray studies, dynam- 
ical analyses, Sunyaev-Zeldovich surveys, and gravitational 
lensing. However, all of these methods require simplifying 



assumptions to be made regarding the shape and/or dy- 
namical state of the cluster in order to derive meaningful 
constraints from available data. Most parametric methods 
typically assume spherical symmetry of the halo and X-ray 
and dynamical estimates additionally assume virialization of 
the cluster. However, examina tion of halo s in C DM struc- 
ture formation simulatio ns (e.g. Bett et al] 1 20071 ) (using the 
Millennium simulation) ; IShaw et al.l 1 200(: ) ) and observed 
galaxy clusters show both of these assumptions to be un- 
physical; simulations show significant triaxiality in cluster- 
scale halos and observed galaxy clusters often exhibit com- 
plex dynamics that suggest r ecent or ongo i ng me rgers and 
disruptio n (e.g. C 10024-H654, ICzoske et al.l l|2002l ): 1E0657- 
558, Clow e et al.l ([2006)). Understanding the impact of these 
physical realities on cluster mass and parameter estimates is 
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crucial for accurate comparisons between measured cluster 
properties and model predictions. 

We focus on the impact of deviations from three as- 
sumptions frequently made in weak lensing analyses of 
galaxy clusters. Gravitational lensing is an appealing tool 
for cluster studies because it is sensitive only to the pro- 
jected mass and thus requires no assumptions be made about 
the dynamical state of the cluster. However, because the 
signal-to-noise ratio is low for weak lensing measurements 
and there is significant degeneracy along the line-of-sight, 
only very simple parametric models can be fit. Even using 
space-based weak lensing data (where the number density 
of background sources useful for a weak lensing analysis is a 
factor of a few higher than that of ground-based data), it is 
often difficult to distinguish between a singular isothermal 
sphere (SIS) with mass density as a function of distance 
from halo centre r give n by p{r) oc an d the univer- 
sal dark matter profile of iNavarro et al.l ([l993) (NFW) with 
p(r) oc in the innermost regions and p(r) oc in the 
outer regions. 

Most cluster profile fits are carried out in the hope of 
either supporting or refuting the universality of the NFW 
profile and thus testing the CDM paradigm. The NFW 
is typically parameterized by an approximate virial mass 
M200 and a concentration parameter, C, and simulations 
predict a strong correlation between the two. For a cluster 
of M = 10^^ M0, C ~ 4. However , sever al authors (e.g. 
[Limousin et al. I l|2007l ): iKneib et all l|2003D : iGavazzi et al] 
( 20031) 1 have recently reported results in the very low prob- 
ability tail of the predicted distribution; notably, in a com- 
bined weak and stron g lensing analysis of Abell 1689, 
iBroadhurst et al.l l|2005l ) report a concentration parameter 
of C = 14 ± 1.5, when C ~ 4 is expected. While one such 
result is not damning, especially given the very complex, 
likely not re l axed, structure of A1689 described recently by 
iLokas et al.l (|2006l ). it is nonetheless of interest to investi- 
gate how possible future discrepancies between observations 
and the predictions of ACDM should be interpreted. 

Crucially, more advanced N-body simulations carried 
out since the ground-breaking work of NFW indicate that 
cluster-scale dark matter halos are expected to be signif- 
icantly triaxial, wit h axis ratios betw ee n minor and major 
axes as small as 0.4 (jShaw et al.l (j2006h ). IOguri et all ()2005l l 
applied a fully triaxial NFW model to the shear map of Abell 
1689 to find that i t is co nsistent with 6% of cluster-scale ha- 
los, and iGavazz 1 (|2005h showed that a triaxial NFW can 
reconcile parameter values derived from observations of the 
cluster MS2137-23 to predictions from N-body simulations. 
Here we carry out a general study of the statistical effects 
of triaxiality on weak lensing analysis. In addition to the 
question of triaxiality, there is also significant scatter in the 
inner and outer slopes of NFW-like profiles in simulations. 
Even the mean val ue of the inner s l ope is contentious, with 
some groups (e.g. iDiemand et al ] liooi); iTasitsiomi et al.l 
l|2004h ') finding it to be steeper than NFW. These effects 
have been explored in various combinations in the con text 
of strong lensing, for example iMeneehetti et all (|2006l l m- 
vestigated the interplay of triaxiality and slope constraints 
in strong lensing analyses. 

We seek to measure the effects of these three alterations 

- triaxiality and variations in the inner and outer slopes 

- to the spherical NFW profile on weak lensing parame- 



ter and mass estimation. Given the difficulty of distinguish- 
ing even entirely different classes of potential dark matter 
profiles using weak lensing data it is not currently feasible 
to meaningfully constrain complex NFW-like models with 
larger numbers of parameters that allow for such variations. 
There are several useful elliptical NFW models obtained by 
perturbing the lensing potential of a spherical NFW model; 
however, these are only realistic for low elli pticities and do 
not account for full triaxial structure (e.g. iGolse fc KneibI 



(|2002h . [Meneghetti. Bartelmann. fc Moscardinil ll2003ll l. We 

therefore estimate the effects of the three expected devia- 
tions from a spherical NFW on the model parameters and 
cluster masses derived fitting the simplest models (SIS, stan- 
dard NFW and the isothermal ellipsoid SIE) to weak lensing 
data. Here we are interested in the properties of individual 
clusters; stacking the lensing signal from a sample of clusters 
gives a much higher signa l-to-noise for the determin ation of 
their average profile (e.g. iMandelbaum et al.l (|2006l )). Also 
note that we consider here the impact of cluster morphol- 
ogy rather than structure along t he line of s ight (for more 
on line of sight structure, see e.g. iHoekstral (120031 )). In the 
next section we introduce the three altered NFW-like pro- 
files and describe their lensing properties as implemented 
in the simulations described in section 3. We present our 
results in section 4, and discuss our findings in section 5. 



2 LENSING BY NFW-LIKE HALOS 
2.1 Weak Lensing Background 

Weak lensing distorts the shapes and number densities of 
background galaxies. The shape and orientation of a back- 
ground galaxy can be described by a complex ellipticity e*, 
with modulus le"! = (1 — b/a)/{l + b/a), where b/a is the 
minomnajor axis ratio, and a phase that is twice the posi- 
tion angle (f), = le^le^**^. The galaxy's shape is distorted 
by the weak lensing reduced shear, g = 7/(1 — ft), where 
7 is the lensing shear and k the convergence, such that the 
ellipticity of the lensed galaxy e becomes 



in the limit of weak defiections. The distributions of ellip- 
ticities for the lensed and unlensed populations are related 
by 



(Pt 



(2) 



assuming a zero-mean unlensed population, the expectation 
values for the lensed ellipticity on a piece of sky is < e > = 
(7 ~ 7. This is the basis for weak lensing analysis in which 
the shapes of images are measured to estimate the shear 
profile generated by an astronomical lens. 

Lensing also changes the number counts of galaxies on 
the sky via competing effects; some faint sources in highly 
magnified regions are made brighter and pushed above the 
flux limit of the observation, but those same regions are 
stretched by the lensing across a larger patch of sky and so 
the number density of sources is reduced. Thus the number 
of sources in the lensed sky n is related to that in the un- 
lensed background no and the slope of the number counts 
of sources at a given flux limit a by n = no^°'~^ , where /i 
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Figure 1. Prolate (a=b=0.4), triaxial (a=0.42, b=0.65), and oblate (a=0.4, b=1.0) clusters of equal mass. 



is the lensing magnification — ( 1 — k)^ — I 7P. A full 
description of these effects is given in ICanizarej () 19821 ). 

We now describe the characteristic behaviour of the con- 
vergence K, and shear 7 of the three analytic NFW-like den- 
sity profiles we study. 

2.2 Triaxial NFW 

A full parameterizatio n for a triaxial NFW halo is given 
bv IJing fc Sutol (|2002l ') (hereon JS02). They generalize the 
spherical NFW profile to obtain a density profile 



p{R) = 



ScPciz) 



(3) 



R/Rsil + R/RsV 

where Sc is the characteristic overdensity of the halo, pc 
the critical density of the Universe at the redshift z of the 
cluster, Rs a scale radius, R a triaxial radius 

= ^ + 5 + ^ ' < ^ ^ = 1)' (4) 
0^ 

and a/c and b/c the minor:major and intermediatemajor 
axis ratios, respectively. In a different choice from JS02 we 
define a triaxial virial radius i?200 such that the mean den- 
sity contained within an ellipsoid of semi-major axis -R200 is 
200pc such that the concentration is 

R200 



C 



Rs 



the characteristic overdensity is 



200 



log(l + C7)-^' 



(5) 



(6) 



the same as for a spherical NFW profile, and the virial mass 
is 

M200 = -^abRiooPc- (7) 

This differs significantly from the parameterization of JS02 
in which the ellipsoidal virial radius is defined in terms of 
an overdensity dependent on the axis ratios. Further, there 
the virial mass is defined in terms of an effective spherical 
virial radius that is a constant fraction of the ellipsoidal 
virial radius, making it independent of the axis ratios of the 
ellipsoid. 

While the parameterization of Jing & Suto has the ap- 
pealing property of giving the virial mass as a function only 



of the virial radius with no need for reference to the axis ra- 
tios of the triaxial halo, it is not ideal for our purposes. Our 
choice of parameters treats the ellipsoid as such all along 
without approximation and gives an effective concentration 
of R2Qo/Rs along each of the three halo axes. Further, our 
choice of an overdensity at collapse independent of axis ratio 
is well motivated by ellipsoidal collapse models that predict 
collapse to st op at the same enclosed density as does spher- 
ical coUapse (jSheth. Mo. fc TormenI (|l999l ')'). Additionally, 
the advantages of deriving a mass from knowledge of the 
virial radius alone are mostly lost in the context of profile 
fitting to observational data, as the axis ratios of the clus- 
ter enter into the characteristic overdensity and thus are 
necessary parameters in any fit of the triaxial model. Most 
crucial to any work is consistency when comparing param- 
eters across different bodies of work; to this end we include 
in Appendix A conversions and comparisons between our 
parameterization and that of JS02, as well as to parameters 
derived by fitting to a spherically averaged density profile, 
similar to those often quoted for N-body simulations. 

Though all our formalism is valid for a general triaxial 
halo, we focus primarily on symmetric prolate and oblate 
halos, the first with two equal-length short axes and the 
second with two equal-length long axes. For ease of distin- 
guishing between these representative halos we introduce an 
axis ratio Q, defined as the ratio between the odd axis and 
the similar axes, so that a prolate halo with a = b — 0.3 and 
c = 1 has an axis ratio Q = 3.33, while an oblate halo with 
a = 0.3 and b = c = 1 has an axis ratio Q = 0.3. Moving 
from small to large axis ratios Q is thus equivalent to begin- 
ning with a very fiat, oblate halo (a "pancake"), stretching it 
along the short axis until it becomes spherical {Q ~ 1), and 
continuing to stretch the same axis past spherical, finally 
ending with a very prolate, "cigar"-like halo. 



2.2.1 Lensing Properties 



prope rties of a triaxial halo 
(|2003| ) (hereon OLS), and 



The full der ivation of the lensing 
is given by lOguri. Lee, fc Sutol , ^ 

we summarize some of that work here|j The triaxial halo 
is projected onto the plane of the sky to find its projected 

^ Although we set c = 1, we keep c as a variable in our notation 
for consistency with OLS. 
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elliptical isodensity contours as a function of the halo's axis 
ratios and orientation angles [6, (j)) with respect to the the 
observer's line-of-sight. The elliptical radius is given by 

e = ^ + ^ (8) 

Ix 1y 

where {X, Y) are physical coordinates on the sky with re- 
spect to the centre of the halo, 



2 



2 



where 
/ = sir 
and 
A = 



2/ 



A + C- ^{A-C)2+B^ 
2/ 



A + C+ y^{A-C)^+B'' 



C 2 
— COS < 



c . 2 
•+ -irsm ( 



+ COS 



C 2 
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H — 5- 7^ sin 



cos 61 sin 2<j!> I - -p: 
b-' 



C 



C 2 

— cos I 



(9) 

(10) 

(11) 

(12) 
(13) 
(14) 



qx 



The axis ratio q of the elliptical contours is then given by 

(15) 

and their orientation angle 'i' on the sky by 
{qx > qy). 



1 -1 B 
* = - tan - — 

2 A- 



C 



(16) 



Here we diverge slightly from OLS's treatment as we are 
interested not in deflection angles but in the lensing shear 
and convergence, both combinations of second derivatives of 
the lensing potential $ (commas indicate differentiation): 



1 



71 

72 



= '^,XY, 
1 



(17) 
(18) 
(19) 



K = -{$,XX+^,Yy). 

These derivatives are calculated as func tions of integrals o f 
the spherical convergence k{(^) (see e.g. iBartelmannI (|l996l ) 
for a full treatment of weak l ensing by a sph eric al NFW 
profile ) following the method of lSchramml (|l990l ) and lKeetonI 
l|200ll '). normalized by a factor of l/\/7 from Equation [TT] 
(see OLS for the derivation of this normalization) 



^,xx 

^YY 
^,XY 

where 



2qX'Ko + qJo, 
2qY^K2 + qJi, 
2qXYKi, 



Kn{X,Y) = 

MX,Y) = 
and 
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1 
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(20) 
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Figure 2. The top (middle) {lower} panels show isoconvergence 
contours for a prolate (spherical) {oblate} halo of M200 = 10^^ 



Mq, C = 4, with axis ratios a 



0.4 {a = b = 1.0) {a 



0.4,6 = 1.0}; the left-hand panel shows the halo oriented with 
the odd axis along the line of sight, the right-hand panel shows 
the halo with odd axis in the plane of the sky. The thick solid lines 
show the limits of the aperture from which weak lensing data is 
taken, and the dashed line shows the Rs ellipse for each projection 
(note that Rs is scaled by the axis ratio in each direction, so that 
when looking at the minor axis of a triaxial halo, the apparent 
scale radius is a times the Rs value of the halo). The lowest 
contour corresponds to k = 0.02 and subsequent contours each 
increase in k by a factor of 2. 



C(^) 



u 

qx 



x^ + 



Y' 



l-il-q')i 



(25) 



Note that our radial variable appears different from 
Keeton's ^ because it is defined in terms of two axis ratios 
qx and qy rather than one q: ( = (,/qx- This reflects a 
dependence on the 3D structure of the cluster; for example, 
extended structure along the line of sight decreases qx and 
thus increases the convergence and shear at a given {X,Y). 

Figure [2] shows the convergence k for a symmetric 
prolate ( "cigar" -shaped), spherical, and oblate ( "pancake" - 
shaped) halo, each of mass M200 = 10^^ Mq. On the left the 
halos are oriented so that the odd axis is along the line of 
sight - the long axis of a prolate halo and the short axis for 
an oblate halo - on the right they are oriented so that the 
odd axis is completely in the plane of the sky. The prolate 
halo oriented with odd axis along the line of sight has the 
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highest convergence due to the large quantity of projected 
mass hidden in it, while the oblate halo in the same con- 
figuration has very low values of n. When viewed with the 
odd axis in the plane of the sky the situation is reversed; the 
prolate halo has low convergence due to a smaller amount of 
mass along the line of sight while the oblate halo has higher 
convergence values. 

Viewing halos in projection is sometimes counterintu- 
itive; for halos of equal triaxial virial mass the 3-D major 
axis of a prolate halo is longer than that of an oblate halo, 
which is in turn longer than the 3-D radius of a spherical 
halo. This is not apparent in projection, where, when the 
odd axis is in the plane of sky, the larger total amount of 
mass in the halo along the line of sight makes the isodensity 
contours of the oblate halo larger than those of its prolate 
counterpart. 



2.3 Varying slopes 

We consider profiles in which only one of the inner or outer 
slopes is allowed to vary from its canonical NFW value. 



2.3.1 Steeper inner slope 

Severa l groups llJirig fc Sutol l|2000l ). iKeeton fc Madaul 
l|200ll ). IWvithe et al.l (120011 )1 have studied a generalized 
NFW with a fixed outer slope p cx r^^ and a varying in- 
ner logarithmic slope — m with density profile 



p{r) 



Scpc 



(r/rs)"" {l + r/rsY 



(26) 



Its lensing properties cannot be computed analytically; the 
convergence is written most simply as 



K{r) 



25cPcrs l-m 

— ^ X 



(3 



) r(y + ^)'"~'(i-yT^)dy 



,(27) 



where x = r/rs and Scr is the lensing critical surface den- 
sity. The she ar is calculated a ccord ing to the relationship 
iMiralda-Escudd (Il99ll )). true for any spheri- 
cally symmetric lens, in which R is the dimensionless mean 
enclosed surface density and is given by 

K(r) = — — X X 



1 



3 — m 
1 



2-Fi [3 — m, 3 — m; 4 — m : —x] + 







iy + xr-^ '-^^ dy\, 

y 



(28) 



where 2F1 is the hypergeometric function. For a spherically 
symmetric lens, the shear is always purely in the tangential 
direction. The characteristic overdensity 5c is given by 



200(3 - m)C" 



32-Fi [3 - m, 3 - m; 4 - m : -C] ' 



(29) 



where C — r2oo s is the concentration pa rameter and r2oo 
is the virial radius (|Keeton fc Madaull200lh FI Note that this 
expression for the overdensity differs by one power of C 
from Equation 2 in Keeton & Madau, correcting a typo 
in the original paper. The virial mass is simply M200 ~ 
|7r200p,rioo. 



2.3.2 Varying outer slope 

The density profile of an NFW-like halo with an outer slope 
left free to vary can be written 

ScPc 



p{r) 



(30) 



(r/r«)(l + r/r,)"-i 

where the exponent n — 1 is chosen so that n — 3 corresponds 
to the canonical NFW logarithmic outer slope of -3. The 
characteristic overdensity is then given by 

_ 200 CHn^2)in~3) 

3 1- [l-h(n-2)C](l + C)2-"' ^ ' 

where C — r2m/i'8', there is no simple expression for the 
convergence or shear so they must be calculated numerically 
by integrating through the mass distribution along the line 
of sight and employing the relationship between R and 7: 



1 



p{\/ r^ + z'^)dz 



'^('^) ~ ~ I r' K,{r')dr' 



10 
R — K. 



(32) 

(33) 
(34) 



2.4 Parametric Fits 

Three simple parametric models are fit to the NFW-like 
models described above. They are the spherical NFW, the 
Singular Isothermal Sphere (SIS) and the Singular Isother- 
mal Ellipsoid (SI E) (a good summary of their lensing be- 
havior is found in lKing fc Schneiderl (|200ll )). They have re- 
spectively 2, 1, and 3 free parameters and are the mod- 
els most often fit to lensing clusters. Each family is char- 
acterized by a set of parameters 11: {C, M200}, {Einstein 
radius 6e}, and {6e, axis ratio q, orientation angle 
respectively. The best-fit parameters are obtained from a 
given lensed catalogue of n.^ galaxies, each with ellipticity 
and position Oj by minimizing th e shear l og-likelihood func- 
tion ( Sch neider, King, fc Erbeiil l|2000l) , iKing fc Schneideij 



(35) 



3 WEAK LENSING SIMULATIONS 

The main body of simulations is carried out for a field 
7.5' in radius, with a background source density no = 



^ Keeton & Madau use a different definition of C— 2 = ^'200/''— 2 
where r_2 is the radius at which the logarithmic slope of the 
density is —2. They give a simple transformation between the 
two definitions, C = (2 — m)C_2- 
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Figure 3. Best-fit concentration C and virial mass M200 for a 
simple NFW model fit to five triaxial lenses with C = 4, M200 = 
10^^ Mq , circularly symmetric on the sky with odd axis oriented 
along the line of sight. Halos with axis ratio Q (ratio between odd 
axis and similar axes) less than one are oblate (Q = 0.3 — ► a = 
0.3), while those with Q greater than one are prolate (Q = 3.3 — > 
a = b = 0.3). 



30/arcminute^, typical of ground-based observations. Al- 
though wide-field imaging of clusters with fields of ~ 0.5 de- 
grees is routinely possible, there are arguments that errors 
due to large-scale structure along the line-of-sight become 
more important as the shear due to the c luster itse l f dimi n- 
ishes with distance from the centre (e.g. iHoekstral l|2003l) ). 
In any case, the general trends will hold for larger fields. 
Poisson noise is accounted for. A catalogue of randomly po- 
sitioned and oriented galaxies with intrinsic shapes e" drawn 
from a Gaussian distribution with dispersion a = 0.2 in the 
modulus |e^| is placed at redshift z = 1. This catalogue of 
background galaxies is lensed through a model lens of choice 
placed at redshift z = 0.18 (the redshift of AbeU 1689), at 
which the width of the field is ~ 1900 kpc/h. Thus our choice 
to place all sources on a sheet at z = 1 is justified by the low 
redshift of our fiducial lens; only for higher redshift lenses 
that are in the heart of the redshift dist ribution is the dis- 
tribu tion of source redshifts important (jSeitz fc Schneid"eil 
1 19971 ). The background galaxies are lensed according to 
Equation[T]and the number counts are reduced as prescribed 
in section 12.11 taking the slope of the source number counts 
in flux to be dlogA''/dlogS' = a = 0.5 (corresponding to a 
slope of 0.2 in magnitude as in lFort. Mellier. fc Dantel-FortI 
l|l997h ). Galaxies located within 1' of the cluster centre are 
removed from the analysis to avoid the strong lensing regime 
at the centre of the cluster (in any case background galax- 
ies near the cluster centre would be mostly obscured by 
cluster members in observations). Throughout we assume 
a concordance cosmology with Q,m = 0.3, Hq — 70 km s~^ 
Mpc~^, h = 0.7, and a cosmological constant S7a = 0.7, 
and a typical massive cluster of triaxial M200 = 10^^ M0 
and C = 4, with corresponding scale radii of Ra ={340.38 
kpc/ft, 461.97 kpc/ft, 626.99 kpc//i} for the spherical, oblate 
(a = 0.4, 6=1), and prolate (a = & = 0.4) cases. At the red- 





Figure 4. The top (bottom) panel shows the mean best-fit NFW 
mass M200 (concentration C) after 500 lensing realizations as 
a function of axis ratio Q for oblate and prolate triaxial halos 
oriented along the line of sight such that the observed density 
contours are circular, for halo concentrations C = 4 and C = 8. 
For small Q the halo is "pancake" shaped and for large Q it is 
"cigar" shaped. Error bars show the approximate Icr dispersion 
within each set of realizations. 



shift of our lens z = 0.18, 1 arcminute corresponds to ~ 127 
kpc//i. 



4 RESULTS 

4.1 Fits to Lensing Data from Triaxial Halos 

To characterize the impact of triaxiality on parameter esti- 
mation over a range of halo shapes we choose ten represen- 
tative triaxial halos to model; five symmetric oblate halos 
with b = c = 1 and five symmetric prolate halos with c = 1, 
a = fe, and a = {0.3, 0.4, 0.6, 0.8, 0.9} for both sets. Each 
halo is studied in two orientations, one in which the odd- 
length axis is oriented along the line of sight (in which case 
an observer would see a circularly symmetric shear pattern 
on the sky), and one in which the odd length axis is oriented 
in the plane of the sky (in which case maximum ellipticity 
is observed) . These will be referred to as the LoS and Plane 
cases respectively. 500 simulations are carried out for every 
halo type and orientation, and NFW and SIE models are fit 
to each of the lensed catalogues. 

We further repeat the simulations using halos of the 
same mass with a higher concentration C = 8 to test how 
sensitive our results are to changes in the underlying cluster 
parameters. Due to the factor of 2 decrease in Rs implied by 
this increase in concentation, we expect the overall quality of 
the fits to decrease as, for some orientations and axis ratios, 
Rs will fall within the 1' annulus removed to avoid the strong 
lensing regime. Throughout, plotted error bars are approx- 
imate la dispersions; though we expect the parameter dis- 
persion within a set of realizations to be non-Gaussian, we 
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use the standard deviation cr^ — (jyTr) ^^iC^^i — 2;)^ as a 
rough indicator of the relative dispersions of fits to different 
halo models. 



^.1.1 Hidden triaxiality: LoS halos 

Figure |3] shows the parameters obtained fitting NFW pro- 
files to 500 le nsing realiza t ions t hrough several of the triax- 
ial halos (see lOguri et alj l|2005l ) for earlier work regarding 
the effects of triaxiality on the M200 — C relationship); in 
Figure U mean best-fit concentrations and virial masses are 
plotted as a function of axis ratio Q for all halos oriented 
along the line of sight. 

We find that LoS oblate clusters are measured to have 
lower masses than true, while prolate clusters are measured 
to have higher masses; in the extreme case of a = 0.3, the 
mass of an oblate halo is underestimated by ~ 40% and 
that of a prolate halo overestimated by ~ 45%. Concentra- 
tions are also significantly affected by hidden triaxiality; for 
extreme oblate halos the measured concentration C is half 
that of the true value while for prolate halos it is double! 
Implied in the concentration and mass fits is a decrease in 
scale radius with increasing prolateness, as expected since 
clusters of equal mass and longer odd axes have shorter cir- 
cular axes; i.e. when one stretches the odd axis, one also 
shrinks the similar axes. Thus for halos oriented with the 
odd axis along the line of sight, the extent of the halo on 
the plane of the sky will be smallest for very prolate halos 
and largest for very oblate ones. 

We find the same type and scale of effects when the 
true concentration is increased to C = 8. We also find sim- 
ilar trends when fitting SIE models, shown in Figure [6l in 
that 9e oc M200 increases with increasing Q: in all cases 
studied a long axis oriented along the line of sight signifi- 
cantly increases mass and concentration estimates, while a 
short axis so oriented reduces them by a similar amount. 



4-1-2 Visible triaxiality: halos in the plane 

Figure [5] plots the mean best-fit concentrations and virial 
masses as a function of axis ratio Q for the studied triaxial 
halos oriented with the odd axis in the plane of sky. Both 
prolate and oblate halos look elliptical on the sky in this 
orientation, as illustrated in Figure [2] In this orientation we 
find that the mass is typically overestimated for oblate ha- 
los (in the most discrepant case of Q = 0.6 by ~ 6%) and 
underestimated for prolate halos (for Q — 3.3 by ~ 45%), 
as expected due to the greater mass hidden in projection in 
the oblate halos. The degree of overestimation of the mass 
decreases slightly for the most extreme oblate axis ratios; 
this can be understood as the effect of increasing elliptic- 
ity on the sky decreasing the amount of mass included in 
the spherical halo fit to the elliptical distribution. The con- 
centration decreases with increasing Q, the opposite of the 
behaviour in the LoS case; however, the deviation form the 
true value is less in this case, with maximum C « 5 and 
minimum C « 3. 

When dealing with significant ellipticity visible on the 
sky it is of course hoped that observations will reveal it and 
better suited models can be fit, perhaps using higher qual- 
ity data. We choose not to use the seemingly obvious choice 





Figure 5. The top (bottom) panel shows the mean best-fit NFW 
mass M200 (concentration C) after 500 realizations as a function 
of axis ratio Q for oblate and prolate triaxial halos oriented with 
the odd axis in the plane of the sky such that the observed density 
contours are elliptical, for halo concentration C = 4. Error bars 
show the la dispersion within each set of realizations. 



of an approximate elliptical NFW for several reasons: such 
models are limited to relatively low axis ratios, cannot give a 
well-defined measure of the 3-D halo mass, and require four 
free parameters, often too many for a weak lensing analysis 
to constrain. Instead, we fit an SIE profile to all ten halos; 
the best-fit SIE Einstein radius and axis ratios q are plotted 
in Figure [G] and A'hoo decrease with increasing Q, re- 
producing the general trend observed in the NFW fits. The 
slight decrease in the best-fitting NFW mass at the low- 
est Q values is not reproduced, further suggesting that the 
decrease observed in the NFW mass is an artifact of the spe- 
cific behaviour of the profile under adverse fitting conditions. 
The SIE generally does a good job of fitting the correct axis 
ratio, plotted in the figure against the actual axis ratio on 
the sky; the mean values of the axis ratios are systemati- 
cally low by ^ 0.05 due to the imposition of a hard wall 
at g = 1 in the minimization routine. Because the offset is 
constant across a large range of axis ratios the SIE emerges 
as an useful tool for detecting visible ellipticity, even when 
the lensing is a result of a non-isothermal distribution. As a 
control we also fit the SIE to LoS halos, all of which have no 
apparent ellipticity on the sky, as shown in the bottom panel 
of Figure |6l We find the same systematic offset of ~ 0.05 for 
the prolate halos; for oblate halos lower values of q are often 
fit with large dispersions due to their low convergence and 
shear values in this orientation. 



4-1-3 Likelihood contours for individual realizations 

Figure [7] plots isolikelihood contours obtained fitting an 
NFW to a single lensing realization through an oblate halo 
with a — 0.4, Q = 2.5, oriented in the plane of the sky and 
along the line of sight, as well as through a spherical halo of 
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Figure 6. The best-fit SIE Einstein radius 9e and visible axis 
ratio qsiE values for C = 4 triaxial halos aligned in the plane 
of sky (along the line of sight) are plotted in the top (bottom) 
panels. For the plane of the sky case, qsiE is plotted as a function 
of the projected triaxial axis ratio q (defined for the triaxial NFW 
in Equation 1151 where q = a, the minor axis ratio, for the plane 
of the sky case). For the LoS cases, for which q = 1, all values 
are plotted as a function of the odd:same axis ratio Q. Error bars 
show the la dispersion within each set of realizations. The C = 8 
halos exhibit similar behaviour, but give systematically higher 
values for Be (the maximum value of 9e is found for the q = 0.3 
prolate LoS case, in which 6e = 46") and for the q = 1 LoS cases 
give qsiE values closer to the true value than do their C = 4 
counterparts. 

the same mass. The contours are tightest for the LoS case; 
in this orientation the halo looks very much like a spherical 
NFW halo with higher mass than true, and so the lensing is 
stronger and the fit more constrained. Conversely, the con- 
tours are larger than for the spherical case when the oblate 
halo is oriented in the plane of the sky, because in this case 
the spherical NFW is a poor fit to the elliptical isodensity 
contours of the lens. The uncertainty of the fit is exacer- 
bated by the lower convergence and shear generated by the 
halo in this configuration. 

4- 1-4 Triaxial Halos Averaged Over Orientation 

In addition to the main body of simulations described above 
additional simulations were undertaken to determine the 
mean distortion caused in parameter and mass estimates for 
triaxial halos across all orientations. For these a higher num- 
ber density n — 120/arcminute^ was used, typical of space- 
based observations. 500 catalogues were lensed through the 
same halo, with the halo at a different random orientation 
for each catalogue. This was done for three oblate halos and 
three prolate halos with a = {0.4, 0.8, 0.9} for each class, 
and for one spherical model as a control. NFW models were 
fit to each catalogue and the best-fit parameters averaged 
across the five hundred orientations. We find that the mean 
parameter values remain very close to the true value, but 



Figure 7. Contours of 70%, 90%, and 99% confidence for indi- 
vidual NFW fits to an M200 = 10^^ Mq , C = 4 prolate halo with 
a = b = 0.4, Q = 2.5, oriented along the line of sight and in the 
plane of the sky. The contours for an NFW fit to a spherical halo 
of the same mass are shown for comparison. Halos with C = 8 
exhibit contours of similar size and shape. 




q q 

Figure 8. The top panels plot the distribution of NFW best-fit 
concentration and virial mass for a prolate halo with a = b = 0.4 
randomly oriented for 500 realizations. The overplotted solid line 
shows the distribution for a spherical halo of the same mass. The 
bottom panels show the distribution of apparent axis ratios on 
the sky of the lensing halos responsible for the two "populations" 
possibly apparent in the prolate halo parameter distributions. 
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that the dispersion around the mean increases with increas- 
ing triaxiality: the dispersion in the mass increases system- 
atically from GM = 0.7 x 10^* in the spherical case 
to rjM = 1-0 X 10^* M0 for the most oblate halo and to 
(J A/ = 2.0 X lO^** Mq for the most prolate halo, and similarly 
the dispersion in the concentration increases from ac = 0.3 
in the spherical case to ctc = 0.8 for the most oblate halo and 
to <Jc = 0.9 for the most prolate halo. This trend is visible 
in Figure |5] where the mass and concentration parameter 
distribution obtained for a spherical halo is overplotted on 
that obtained from an a = 6 = 0.4, Q = 2.5 prolate halo of 
the same mass. 

The distributions for the prolate halo are significantly 
more non-Gaussian than their spherical counterparts; even 
more, the distribution of C values appears somewhat bi- 
modal, with most halos falling in a Gaussian-like distribu- 
tion around the true value C = A but a significant number in 
a subpopulation centered at C « 5.8. Plotting the axis ratios 
on the sky q (as calculated in Equation I15|) of the member 
halos of these two apparent "populations," as is done in the 
bottom panels of Figure [S] shows that the sub-group on the 
right contains most of the g ~ 1, near LoS orientated lenses. 
It seems then that halos with significant triaxiality visible 
in the plane of the sky are effectively radially averaged by 
fitting with a spherical model, and thus the parameters fit 
are evenly distributed near the true value, with a slight skew 
to low values. However, halos that appear almost spherical 
on the sky are consistently fit as more massive, more con- 
centrated halos in this case of an underlying prolate lens. 
Oblate halos show the opposite behaviour, leading to con- 
sistent mass and concentration underestimates when LoS 
oriented. This behaviour is illustrated for the a — 0.4 pro- 
late halo in Figure [S] which plots parameters obtained from 
a single lensing realization under very low noise conditions 
in which the number density of sources no is very high and 
the ellipticity dispersion a is very low, as well as a scaled 
likelihood (0 best, 1 worst) as a function of orientation angle 
9. As expected, the likelihood decreases significantly as the 
halo moves from LoS to Plane orientation due to the increas- 
ing poorness of fit of a spherical model to the increasingly 
elliptical profiles. 



4-1.5 Behavior of parameters as a function of Q 

Catalogues simulated with very low ellipticity dispersion a 
and high number density no minimize the effects of noise, 
generating, for Gaussian-like error distributions, best-fit pa- 
rameters that represent the mean best-fit parameters that 
would be obtaine d after running many reali z ations for real 
noise conditions l|Schneider. King, fc ErbenI (|200G| )). Cata- 
logues were lensed under such reduced-noise conditions for 
many values of a in both oblate and prolate triaxial mod- 
els aligned with the odd-axis along the line of sight, to un- 
derstand the "worst-case" behaviour of parameter estimates 
across a range of axis ratios. Figure [TOl shows that the mass 
and concentration estimates increase with increasing Q\ the 
mass function fiattens towards extreme prolate axis ratios, 
while the concentration function remains steeper. This sug- 
gests that when very high concentrations are measured for a 
given mass, a halo may be in this regime of significant elon- 
gation hidden along the line of sight. A scaled likelihood 
value for each halo is also plotted (0 corresponds to the fit 
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Figure 9. Best-fit NEW concentration, virial mass, and scale 
radius and the scaled likelihood of those fits for a prolate halo 
{a = b = 0.4) rotated through 90° on the sky. The dotted lines 
show the true parameter values for the lensing halo. The best-fits, 
highest concentrations, and highest masses are obtained when the 
long axis of the halo is oriented along the line of sight. 
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Figure 10. Best-fit NFW concentration, virial mass, and scale 
radius Ts and the scaled likelihood of those fits for halos aligned 
along the line of sight as a function of axis ratio Q. Boxes indicate 
the location of the spherical model. 



to a true spherical model, 1 is the worst fit in our sample); 
oblate halos are fit quite poorly compared to their prolate 
counterparts due to their lower convergence and shear in the 
LoS orientation. 
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Figure 11. The top (bottom) panel plots the mean best-fit 
NFW mass M200 (concentration C) after 500 lensing realizations 
through halos with variable inner slope, as a function of the inner 
slope m, for halo concentrations C = 4 and C = 8. Error bars 
show the 1(7 dispersion within each set of realizations. 



Figure 13. The top (bottom) panel plots the mean best-fit 
NFW mass M200 (concentration C) after 500 lensing realizations 
through halos with variable outer slope, as a function of outer 
slope parameter n, for halo concentrations C = 4 and C = 8. 
Error bars show the Icr dispersion within each set of realizations. 



mean concentration Cfu ~ 15, and a halo with C = 8 is 
fit with mean concentration Cfu ~ 24). Thus assuming the 
paradigmatic value for the inner slope when in fact there 
is uncertainty about its mean value and significant scatter 
about that mean can have serious implications for the dis- 
tribution of expected observed concentrations. The mass is 
underestimated for steeper-sloped halos, and so therefore is 
0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 the scale radius. The SIS models exhibits the opposite trend, 

^ preferring larger 6e values and masses for lenses with steeper 

inner slopes. This likely occurs because the NFW-like halo 
Figure 12. The mean best-fit SIS Einstein radius Oe after 500 ig becoming more like an SIS as its inner slope approaches 

lensing realizations through halos with variable inner slope, as a isothermal, which has a naturally higher Einstein radius at 

function of m, for concentrations C = 4 and C = 8. Error bars 

' a given mass, 

show the la dispersion within each set of realizations. 
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4.2 Fits to Lensing Data from Halos with Varied 
Slopes 

4-2.1 A steeper inner slope 

Six NFW-like halos with varied inner slope are studied with 
values of m = {0.8, 1.0, 1.2, 1.4, 1.6, 1.8}. Though there are 
no predictions calling for inner slope values of less than one, 
m = 0.8 is included to account for scatter about the mean 
and to allow observation of trends across the canonical NFW 
value m = 1. 500 catalogues of background galaxies are 
lensed and NFW and SIS models are fit to the simulated 
lensed catalogues. We further repeat the simulations using 
halos of the same mass with a higher concentration C = 8 
to test how our results differ with changes in the underlying 
cluster parameters. 

The best-fit NFW and SIS parameters are shown in Fig- 
ures [TT] and [T51 We find that NFW-like halos with steeper 
inner slopes give rise to very high estimated concentrations 
(for m = 1.8, a cluster with a true C = 4 is fit with a 



4-2.2 Scatter in the outer slope 

Seven NFWs with varying outer slopes are simulated, each 
with C = 4 and M200 = 10^^ M©, with values of n = 
{2.4,2.6,2.8,3.0,3.2,3.4,3.6} to understand the effects of 
scatter about the canonical value on lensing analyses. 500 
catalogues are generated and lensed as described above for 
each halo model and NFW and SIS models are fit to each. 
Again, we repeat the simulations with a concentration of 
C = 8. 

The best-fit NFW and SIS parameters are shown in Fig- 
ures [TJ] and [TJ] The NFW concentration is overestimated for 
steeper outer slopes and underestimated for shallower slopes; 
the mass follows the concentration, unlike in the case of a 
steeper inner slope. The mass is a nearly linear function of 
steepness; an increase in steepness leads to a mass increase 
comparable to the decrease in M200 caused by a slope de- 
crease of the same magnitude. The behaviour of the concen- 
tration, however, is not linear, with increases in steepness 
producing larger changes in C than comparable decreases. 
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Figure 14. The mean best-fit SIS Einstein radius 9e after 500 
lensing realizations through halos with variable outer slope, as 
a function of outer slope parameter n, for halo concentrations 
C = 4 and C = 8. Error bars show the la dispersion within each 
set of realizations. 



This means that even slopes normally distributed around 
the paradigmatic value n — 3 will lead on average to mea- 
surements of slightly higher concentrations. The SIS fits dis- 
play the same mass trend, with mass and 9e increasing with 
increasing steepness. 
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Figure 15. The distribution of best-fit NEW concentration and 
virial mass for 10,000 lensing realizations through 10^^ Mq halos 
with shapes drawn from the simulation-predicted triaxiality dis- 
tribution and randomly oriented. The solid line plots the same 
population, this time including only realizations for efficiently 
lensing halos with < g > ^ 0.014, the average shear for a 10^^ 
Mq spherical halo. 



5 DISCUSSION &: CONCLUSIONS 

We focus our further discussions on the impacts of triaxi- 
ality, the most potentially insidious of the three deviations 
from NFW we study, as it cannot be fully investigated and 
marginalized by any lensing technique due to degeneracies 
along the line of sight. 

To begin, it is crucial to know how much triaxiality is 
ex pected in the clus ter-sized halos. Recent N-body studies 
bv lshaw et all l|2006l ) have resulted in detailed predictions of 
the axis ratio distribution for dark matter halos, finding that 
halos tend to be more prolate than oblate, with the distribu- 
tion for b/c peaking near 0.8 and that for a/c peaking near 
0.65 with a tail down to 0.4. While the inclusion of baryon in- 
teractions might reduce the levels of predicted triaxiality (H. 
Hoekstra, private communication), Shaw's predictions are a 
very useful limiting case in which to apply our results. The 
top panels of Figure [TJ] plot the parameter distributions ob- 
tained from performing 10,000 lensing realizations through 
a 10^^ M0 , C — 4: halo, at each realization choosing axis ra- 
tios from the predicted distributions and randomly orienting 
the halo. The means are equal to the true values, and the 
distributions are close to Gaussian, but with very long tails 
toward high values. Though those tails are indeed extended 
by triaxiality, indicated by small increases in the dispersion 
of the distributions from those for distributions obtained by 
lensing through spherical halos, the total effect of triaxiality 
on a complete population of clusters is very small. 



5.1 Lensing Efficiency Effects 

Our finding that the mean best-fit C and AI20Q approach 
the true value for even very triaxial halos when averaged 
over orientation initially suggests that while triaxiality may 
increase the errors on parameter estimates for individual 
clusters, it should not be a significant issue in constraining 
the mass function when tabulated over a large sample. How- 
ever, this conclusion neglects to take into account the effects 
of orientation on lensing efficiency. 
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Figure 16. Average reduced shear g within the annulus located 
between 1' - 7.5' of the halo center, the annulus from which the 
weak lensing data is drawn in the simulations, plotted as a func- 
tion of axis ratio Q for LoS and Plane halos. 



Certain orientations of clusters will produce higher lens- 
ing efficiencies than others, making them more likely to be 
included in lensing-selected cluster surveys. Figure [161 plots 
the average reduced shear over the annulus from which the 
weak lensing data is taken as a measure of lensing efficiency. 

Halos with significant mass along the line of sight are 
the best lenses; prolate halos with high Q (low a and b) 
oriented along the line of sight have almost twice the average 
shear than does a spherical halo of the same mass. Oblate 
halos with very low Q (low a and b) oriented in the plane 
of sky are also effective lenses, which can be understood 
by recalling that the long axis visible on the sky is also 
the length of the axis hidden in projection, giving rise to 
significant convergence and shear. 

Halos that are the most efficient lenses will be more 
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likely to be studied as lenses; thus clusters selected for their 
lensing properties are likely to reside in the high-triaxiality, 
low probability tails of the halo distribution. While triax- 
ial halos averaged over orientation in simulations give mean 
parameter values equal to the true values, this type of sim- 
ulation does not recreate observational conditions in which 
prolate halos are much more likely to be observed as lenses 
when oriented near the line of sight, and oblate halos when 
oriented in the plane of the sky. Furthermore, since very 
prolate and oblate halos (aligned advantageously) are much 
better lenses than spherical halos, observed halos, especially 
those with strong lensing effects such as Abell 1689, are 
likely to be represent the high triaxiality end of the halo 
distribution. 

Understand ing this behavi o r is i mportant for survey 
design (see e.g. IWittman et al.1 l|2006l ) for more on shear- 
selected surveys); any criteria involving lensing efficiency 
will skew both mass and concentration estimates high! As an 
example, the bottom panels of Figure [TJ] shows distributions 
of NFW parameters for the same population of 10,000 lens- 
ing realizations shown above, this time keeping only param- 
eters from realizations in which the lensing halos have more 
than a cutoff level of lensing efficiency, set to be the lensing 
efficiency of a spherical halo of 10^^ Mq, < g >= 0.014. 
Though even after this cut the effects of triaxiality are very 
small, the mean parameter values do indeed shift slightly 
towards higher mass and concentration. 



5.2 Model Discrimination 

A potential impact of the deviations from NFW that we 
study here that we have not yet addressed is that of changing 
the frequency with which one model family is incorrectly 
preferenced over another. We find that varying the outer 
slope has no significant impact, and that triaxiality has only 
a very small impact, with very triaxial halos of a = 0.4 being 
misidentified as SIS halos in 4% more realizations than are 
spherical halos. However, changes in the inner slope prove 
to be very significant; for an inner slope of m = 1.6 NFWs 
are misidentified as SIS halos in 83% more realizations than 
are canonical NFW halos. Interestingly, an NFW with m — 
1.8 fares significantly better, with misidentification in only 
29% more realizations. It is not surprising that as the inner 
slope of the NFW model approaches the isothermal slope 
it would be confused more often with an SIS, and that fact 
that it fares worst when m = 1.6 suggests that the SIS 
fitting routine is "averaging" the slope in some sense over the 
range in which the data is taken. It is thus clear that weak 
lensing analyses are quite sensitive to the inner slope, even 
though the data are taken from outside the inner regions of 
the cluster, due to the non-local nature of the gravitational 
potential. 

5.3 The Meaning of Concentration 

Our choice t o use the traditional C = r2oo/ra, as did 
IWvithe et al] (|200ll '). rather than C-2 of iKeeton fc Madaul 
l|200ll ) in our treatment of an NFW with a free inner slope 
derives from two factors. The first is that C is the more 
commonly used parameter and t hus more useful fo r com - 
parisons; the second is that while iKeeton fc Madaul (|200ll ) 



argue that the point at which the logarithmic slope is —2 is a 
more physical quantity than the scale radius , we feel that 
the break radius at which the slope of the NFW is the av- 
erage of the asymptotic inner and outer slopes is an equally 
valid physical scale to use. In any case, the numerical val- 
ues of concentration are most important in comparison to 
those from simulations and between different members of 
the cluster population, using the same definition. 



5.4 Summary and Outlook 

In this paper we concentrate on constraints from weak lens- 
ing; much work remains to be done, marginalising over 
the effects of t r iaxial ity for cluster data sets, pioneered 
by lOguri et al] (120051 ) for AbeU 1689. In the past cou- 
ple of years, a great deal of progress has been made in 
techniques combining weak and stro ng gravitational lensing 
constraints o n clu ster profiles (e.g. ICacciato et al.l l|2006l ): 
iDiego et al.l (|2007t )). In addition, gravitational flexion has 
very recently been harnassed by [Leonard et al.l (|2007l ) to 
probe substructure in clusters, and is also a promising 
means of constraining the mass distribution in the regime 
where standard weak lensing approximations break down. 
Besides the various regimes of gravitational lensing, large 
spectroscopic surveys have the potential to provide infor- 
mation essential to constructing models of the 3-D distribu- 
tion of cluster mass and its dynamical history. So far this 
has only been undertaken fo r a very small number of clus- 
ters (e.g. lCzoske et al.l l\200t l). but the advent of wide-field 
spectroscopy makes this feasible. Such spectroscopic stud- 
ies combined with deep space-based observations (where the 
number density of galaxies useful for shear measurements is 
several times that of ground-based data) and SZ and X-ray 
studies will allow for more detailed modeling of the mass 
distributions of clusters. 

We investigate the effects of triaxiality and variations 
in the slopes of the NFW profile on weak lensing parameter 
and mass estimates and find that for individual clusters the 
effects can be significant. Triaxiality causes the mass and 
concentration to be over or under estimated by ~ 50% and 
a factor of 2, respectively; steepening the inner slope leads 
to up to 10% underestimation of the mass and overestima- 
tion of the concentration by a factor of up to 3; scatter in 
the outer slope causes up to 10% and 5% errors in mass and 
concentration estimates, respectively. Averaged over orien- 
tation and over predicted axis ratio distributions the errors 
induced by triaxiality are significantly reduced, which bodes 
well for the use of cluster masses as a cosmological probe; 
however, when lensing efficiency is taken into account they 
become more important. Thus the selection effects must be 
accounted for in survey design, for example when shear se- 
lected samples are employed or in future follow-up of clusters 
detected using the S-Z effect. Crucially, the effects of triax- 
iality along the line of sight and variations in halo slopes 
may lessen tensions between the unusual galaxy cluster con- 
centration/mass estimates reported in recent work and the 
predictions of ACDM. 
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APPENDIX A: TRIAXIAL PARAMETER 
CONVERSIONS 

We here compare o ur triaxial NFW parameters to those of 
IJing & Sutol Jigoi) as well as parameters fit to spherically 
averaged 3D ellipsoidal density profiles, as is often done in 
fits to clusters in N-body simulations. 

In Figure lAll the JS02 concentration parameter Cjs 
and virial mass Mjs are plotted against C for a clus- 
ter of M200 = 10^^ M0 for two prolate clusters and a 
spherical NFW. As expected, there is close agreement be- 
tween the parameters when the clusters are nearly spher- 
ical, but there are significant differences at more extreme 
axis ratios. To transform between parameter systems we re- 
quire that a halo described in our coordinate system by 
n = {C,Rs,R20o,M2oo,a,h} give the same density profile 
as a model described by H' = {Cjs, Rq, Re, Mjs,a,h} in 
JS02's parameterization. Cjs is defined as 



Cjs '■■ 



Re 

R^ 



(Al) 



where Re is defined such that the mean density within an 
ellipsoid of semi-major axis Re is 
0.75 

P — 5Avir ( ^ 



(A2) 



where A„ir is the overdensity at virialization of the halo. 
The characteristic overdensity of the halo. See is then given 
by 



5A„ir I Jji, 



ln(l + Cjs) 



Cj 



(A3) 



14-Cjs 



The second term is the same as for our parameterization, 
while the prefactor depends on the axis ratios unlike in our 
model. Since Tio = Rs is the same in both systems, setting 
the density profiles equal simply requires that 5e = See for a 
given set of parameters. Thus, converting from our param- 
eterization to that of JS02 requires the numerical solution 
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Figure Al. The JS02 concentration parameter Cjs and virial mass Mjs plotted as a function of C for lialos of Af200 = 10^^ Mq. 
Three examples are shown; one a very prolate halo with a = b = 0.4, one a slightly prolate halo with a = b = 0.8, and a spherical halo. 
In the JS02 parameterization the halo is characterized as less concentrated and more massive for all triaxial models. 
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Figure A2. The spherical concentration parameter Cgphere and virial mass M^pf^^^.^ as a function of triaxial axis ratio Q for clusters of 
.'^^200 = 10^^ Mq and C = 4. The JS02 concentration and mass are also shown for comparison. 



of 



ln(l + Cjs) 



Cj 



1 + Cjs 



200 



(A4) 



Re is then simply calculated, Re = CjsRo, and Mjs is 
given by the fitting formulae derived from JS02's numerical 
simulations: 



Mjs = 



4:TTAviM{z)pe{z) 



f— v 

Vo.45y 



(A5) 



Figure lA2l shows the best-fit spherical concentration pa- 
rameter C sphere and virial mass M^phere as functions of axis 
ratio Q, defined in Section [4. 1.1 1 Parameters are calculated 



by averaging the density of the ellipsoidal cluster over spher- 
ical shells and fitting a spherical NFW density profile to the 
averaged values. This is how fitting is often done to N-body 
simulations, in which parameter distributions are tabulated 
assuming spherically symmetric halos due to the significant 
degeneracies encounte red in fitting triaxial density models 
(|Diemand et al.1 l|2004l )'). The spherical concentration is very 
similar to our choice of C at most axis ratios, but the spher- 
ical mass is significantly less than our choice of ellipsoidal 
M200 for very large and small Q. 

This paper has been typeset from a T^^/ WTfi]K. file prepared 
by the author. 



